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ABSTRACT 

High-resolution circular spectropolarimetric observations, obtained with ES- 
PaDOnS in the context of the BinaMIcS Large Program, have revealed a magnetic 
held in the B3V secondary component of the SB2 binary system e Lupi (B2/B3). As 
the B2V primary is already known to be magnetic, this is the hrst detection of a 
magnetic held in both components of an early-type binary system. The longitudinal 
magnetic held of the primary is ~ —200 G; that of the secondary ~ -1-100 G. Observa¬ 
tions can be approximately reproduced by a model assuming the magnetic axes of the 
two stars are anti-aligned, and roughly parallel to their respective rotation axes. Esti¬ 
mated magnetospheric radii indicate a high probability that their magnetospheres are 
interacting. As many of the arguments for the different proposed formation scenarios 
of fossil magnetic helds rely upon evidence drawn from investigations of close bina¬ 
ries, in particular the rarity of magnetic ABO stars in close binaries and the previous 
absence of any known close binary with two magnetic, massive stars, this discovery 
may be an important new constraint on the origin of fossil magnetic helds. 
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1 INTRODUCTION 

A principal result of the Magnetism in Massive Stars 
(MiMeS) snrvey has been the demonstration that a sig- 
nihcant snbset (~10%) of the early-type stellar popula¬ 
tion host strong (0.1-10 kG), stable magnetic helds (Wade 
et ah, submitted). In contrast, magnetic helds are essen¬ 
tially ubiquitous amongst cool stars. The magnetic helds 
of magnetic, massive stars ( MMSs) are cur r ently explained 
as so- called fossil helds (e.g. ICowlind Il945l : iMestel fc MossI 
I 1977 I : iBraithwaite fc SDruitll2004h : magnetic hux preserved 
from a previous era in a star’s life. Magnetohydrody¬ 
namic (MHD) simulations indicate that magnetic helds 
within a radiative envelope can relax into stable conhg- 
urations which then dissipate due to Ohmic decay and 
other slow mechanisms on stellar ev o lution ary timescales 
llMestel fc Takh^ Il972l : iBraithwaitel I 2 OO 9 II . This frame¬ 
work has been highly successful in explaining the ob- 
served characterist ics of the magnetic helds of hot stars 
llNeiner et al jlmi), predicting magnetic topologies domi¬ 


nated by strong, global dipoles, wi th large angles b e tween 
magnetic and rotational axes (e.g. iDuez et al.l[201(]|. I 2 OIII: 
lEmeriau fc Mathid l2015h . These characteristics are essen¬ 
tially identical to those long kno wn for the magnetic Ap stars 
(e.g. iBorra fc Landst reed 198dl and the He-weak and He- 
stron g Bp stars (e.g. iBorra fc Landstreetlll979l : iBorra et al.l 
Il983ll . suggesting a common mechanism behind the mag¬ 
netic helds of all stars with radiative envelopes, i.e. from 1.5 
to 50 Mq llWadell2015l ~l. 


Notwithstanding these successes, the origin of fossil 
helds remains a subject of speculation, indicating a fun¬ 
damental dehciency in our understanding of a basic phys¬ 
ical property of all stars more massive than 1.5 Mq. One 
of the principal goals of the Binarity and Magnetic Interac¬ 
tions in various classes of Stars (BinaMIcS) project has been 
to investigate the o rigin of fossil helds amongst the MMSs 
llAlecian et al.ll2015l L by examining in particular the prop¬ 
erties and prevalence of binary systems containing MMSs. 

One of the BinaMIcS targets is e Lupi (HD 136504, 
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B2/B3), an SB2 system comprised of 2 early B-type stars 
with an orbital period of ~4.6 d. Tentative e vidence for a 
magne tic field in the system was presented bv iHnbrig et alJ 
(l2009tl . based upon FORSl data. The first reliable detection 
was confirmed with ESPaDOnS data bv IShultz et al.l ll2012ll . 
Here we report for the first time the detection of a magnetic 
held in both components of e Lupi, making it the hrst known 
close binary with two magnetic, massive stars. 


2 OBSERVATIONS 


Our dataset is comprised of 51 Canada-France-Hawaii Tele¬ 
scope (CFHT)/ESPaDOnS spectropolarimetric sequences. 
ESPaDOnS is a high-resolution (A/AA ~ 65,000) spec- 
tropolarimeter with spectral coverage from 369.3 nm to 1048 
nm across 40 overlapping spectral orders. Each sequence is 
comprised of 4 polarized subexposures, which are combined 
to yield unpolarized intensity (Stokes I), circularly polarized 
hux (Stokes V), and diagnostic null {N) spectra. On each 
night a minimum of 4, and a maximum of 11, sequences 
were obtained. Sequences obtained on the same night were 
co-added in order to increase the SNR, yielding 10 inde¬ 
pendent measurements. The hrst two were acquired by the 
MiMeS Large Program (LP), the most recent by the BinaM- 
IcS LP, and the remainder in the context of PI program 
CFHT14AC010 (PI: M. Shultz). The log of observations is 
provided in Table [T] 

To further increase the SN R, we used the leas t-squares 
deconvolution procedure fLSD:lD onati et al.||l997lf. in par¬ 
ticular the iLSD package dKochukhov et ^ 2010h . The line 
list used to construct the LSD mask was obtained from 
the Vienna Atomic Line Database VALD3 llPis kunov e t al.l 
I 1995 I : iRvabchikova et al]ll997l : iKupka et al.lll999l . l2000ll . us- 


ing an ‘extract stellar’ request. Since both components have 
similar temperatures, we used a single mask that we ad¬ 
justed to best match the spectrum. We began with a 20 kK 
solar metallicity mask, cleaned to include only metallic lines 
unblended with H, He, or telluric lines, with 234 spectral 
lines remaining in the mask after cleaning. The remaining 
lines were ‘tweaked’ by empirically adjusting the line depths 
to match the observed line strengths. LSD profiles were then 
extracted using 3.6 kms“^ velocity pixels, with normaliza¬ 
tion values of the wavelength and Lande factor of Ao = 500 
nm and go = 1-2. Noise was suppressed using a Tik honov 
regularization factor of 0.2 dKochukhov et ahl I2OI0II . The 
LSD profile extracted from the BinaMIcS observation (the 
highest SNR measurement, obtained on HJD 2457122) is 
shown in Fig. [T] 

To quantify the presence of a magnetic signature, we 
applied two diagnostic techniques to the LSD prohles. The 
hrst was to calculate the False Alarm Probabi lity (FAP) of 
the s ignal in Stokes V inside the line prohle dPonati et al.l 
Hii^). Detections are considered dehnite (D D) if FAP < 
10“®, and non-detections (ND) if FAP > 10“® dPonati et al.l 
[T^. Integration ranges were set individually for each LSD 
prohle, depending on the radial velocities and line-prohle 
widths of the components, as illustrated in Fig.[T] Detection 
hags are provided in Table [T] The most recent, high-SNR 
measurement yields a DD in the secondary component. 

The longitudinal magnetic held (Bz) was measured 
from the LSD prohles by taking the hrst-order moment of 



Figure 1. LSD profile from the night of 2015/04/09. Dotted lines 
indicate the integration range for the primary, dot-dashed lines 
for secondary. Both stellar components are separated in Stokes 
I (bottom); N is consistent with noise (middle); in Stokes V, a 
definite detection is registered in both line profiles (top). 


the Stokes V prohle (lMathvElll989f ). The same integration 
ranges were used for measuring (Bj,) as for evaluation of the 
FAPs. (Bz) values for combined line prohles are provided 
at all orbital phases in Table [T] and for the primary and 
secondary components individually at phases at which the 
components’ lines are separated. We also measured (Az) us¬ 
ing the N prohles; it is consistent with noise in both the 
blended prohles and the individual prohles at all phases. 

Prior to the clear detection of the magnetic signature 
of the secondary in April 2015, analysis of data collected 
between 2011 and 2014 suggested that the secondary com¬ 
ponent might also possess a magnetic held. At phases in 
which the components are clearly separated, the secondary 
consistently registers (Bz) ~ -1-100 G, while the magnetic 
held of the primary is always ~ —200 G. When the compo¬ 
nents are blended, (Bz) ~ —100 G, as might be expected 
if two Zeeman signature of opposite polarity are partially 
cancelling. All 7 FO RSl measurements are also ~ —100 G 
jHubrig et al.l 1^1 1^ : due to the low spectral resolution of 
this instrument, the components are likely blended at all 
orbital phases. 


3 DISCUSSION 

e Lupi’s orbital eleme n ts we re most recently determined 
by lUvtterhoeven et al.l |(200l) (hereafter U05): e = 0.277, 
uj — 18°, Borb ~ 4.56 d,i = 21°. The LSD prohles are 
shown phased with the U05 ephemeris in Fig. U05 found 
usinf = 42 kms“^ for the primary and 37 kms“^ for the 
secondary. In the left panels of Fig.[2]we show model Stokes 
7 LSD prohles. The basic disk integration model used in 
synth esizing the LSD prohles is described bv IPetit fc Wad3 
ll2012l) . The model approximates the LSD prohle as an in¬ 
dividual spectral line, rather than extracting an LSD pro¬ 
hle from a synthet ic spectrum; while ther e are limitations 
to this approach dKochukhov et ahl I2OIOI) , it is sufficient 
for a first-order model of the observed LSD profiles, and 
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Table 1. Log of observations and magnetic measurements. SNR is per 1.8 kms“^ spectral pixel. DF refers to the detection flag, as 
discussed in the text. Columns 7—9 (‘Combined’) provide magnetic measurements performed including both stellar line profiles; at phases 
during which the profiles are blended, ony combined measurements are given. Remaining columns give magnetic measurements for the 
individual stellar components, at phases at which the components’ line profiles are well-separated. 


HJD - 

2450000 

UT Date 

*5^orb 

SNR 

Combined 

(Bz) {N,) 

(G) (G) 

DF 

(Bz> 

(G) 

Primary 

{Nz) 

(G) 

DF 

Secondary 
(Bz) (iVz> 

(G) (G) 

DF 

7122.0250 

2015-04-09 

0.996 

4904 

-102T30 

-43±30 

DD 

-183T12 

-16T12 

DD 

104T14 

13±14 

DD 

6824.8751 

2014-06-16 

0.827 

1642 

-127T42 

4±42 

DD 

-240±25 

-28±25 

DD 

78T29 

8±29 

ND 

6822.8914 

2014-06-14 

0.392 

1343 

-119T60 

-39±60 

DD 

-228±37 

-59T37 

DD 

104±33 

-51±33 

ND 

6821.8707 

2014-06-13 

0.168 

1212 

-102T28 

-7±28 

DD 

- 

- 

- 

- 

- 

- 

6819.8677 

2014-06-11 

0.729 

2239 

-115±13 

19T13 

DD 

- 

- 

- 

- 

- 

- 

6816.8524 

2014-06-08 

0.067 

1779 

-196±45 

-65±45 

DD 

-220±24 

-16T24 

DD 

83±29 

2±29 

ND 

6760.9572 

2014-04-13 

0.809 

1326 

-103±37 

-8±37 

DD 

- 

- 

- 

- 

- 

- 

6756.9622 

2014-04-09 

0.933 

1261 

-94T95 

-100T95 

DD 

-175±44 

-33T44 

DD 

147±46 

-57±46 

ND 

5727.8117 

2011-06-15 

0.227 

1461 

-127±33 

19±33 

DD 

- 

- 

- 

- 

- 

- 

5634.1538 

2011-03-13 

0.687 

1325 

-124T22 

-17±22 

DD 

- 

- 

- 

- 

- 

- 
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Figure 2. Observed (solid black lines) and synthetic (dashed 
red lines) LSD profiles arranged in order of orbital phase. The 
BinaMIcS observation is at phase 0.996. 

much cheaper computationally. Due to the inclusion of an 
ad hoc turbulent broa dening, chosen to op timize the fit to 
the observed profiles JPetit fc Wadell2012l) . we use slightly 
lower usini values than those given by U05: u sin ip = 37 
kms“^and u sin is = 27 kms“^. The synthetic LSD profiles 
corresponding to the individual stars were centred at radial 
velocities computed using the U05 orbital elements, and then 



- 0.2 - 0.1 0.0 0.1 0.2 
X (AU) 


Figure 3. Schematic of the e Lupi system’s orbit, in the inertial 
frame of the primary. The primary is indicated by the solid (red) 
circle in the centre, the secondary by the solid (blue) circle at 
periastron. The orbit of the secondary is indicated by the solid 
(blue) line. The dashed (red) line indicates the primary’s R\. The 
dot-dashed (blue) lines indicate the inner and outer extent of the 
secondary’s Rx during an orbital cycle. Orbits and stellar radii 
are shown to scale. 

added using a flux ratio of 70%/30% primary/secondary, 
matching the observed line strength ratio. The combined 
profiles were then normalized to reproduce the equivalent 
widths of the observed LSD profiles. 

The right panels of Fig. [2] show a comparison between 
observed and synthetic Stokes V profiles. While the rota¬ 
tional periods, inclinations i of the rotational axes, obliq¬ 
uity angles j3 between the rotational and magnetic axes, 
and dipolar field strengths Ba cannot yet be determined 
with certainty, some inferences can be made from the ex¬ 
isting dataset. At phases during which the profiles are well- 
separated, (Bz) is constant within the error bars, while at 
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blended phases, (-Bz) is always about —100 G. This suggests 
that {Bz) of both components is roughly constant, implying 
that either i, /3, or both, are relatively small (considering the 
known v sin is, the other possibility - that the periods are 
much longer than the timescale over which the data were 
acquired - is not considered). We therefore set i to the or¬ 
bital inclination of 21 ° for both stars, and fixed j3 = 0° 
such that all variation in Stokes V is due to orbital modula¬ 
tion, consisting only of RV shifts applied to constant Stokes 
V profiles. We adopted Ba.p = 900 G and Bd,s = 600 G, 
as being approximately the minimum strengths compatible 
with the observed (Bz). 

As is clear from Fig. (2] this model reproduces the gross 
features of Stokes 7 and V at all phases, including those dur¬ 
ing which the line profiles are blended. In many cases, small- 
scale features are present in the observed Stokes V pro- 
files. These are like ly a result of non-radial pulsation (U05, 
iTelting et al.l l2006l l. as indicated by comparisons between 
the observed Stokes I and V profiles: strong asymmetries in 
I are typically accompanied by corresponding asymmetries 
in V, as is seen most clearly at phase 0.392. 

Determining the rotational periods and magnetic prop¬ 
erties of these stars will require a larger, high-SNR spec- 
tropolarimetric dataset. Modelling the nonradial pulsations 
will be necessary, as these effects may well dominate the 
variability of Stokes V for each star; as these frequencies 
and modes are unidentified (U05), large, high-quality pho¬ 
tometric and spectroscopic datasets will also be needed. 

We evaluate the strength of magnetic wind confinement 
vi a the magnetic wind confi nement parameter 77 * (Eqn. 7 
oflud-D oula, &: Owockil [2002l i . using the mass-loss recipe of 
IVink et al.l I 200l| i. Taking 7?,,p = 4 . 9 ±O. 8770 , M*,p = 8.7± 
I.OMq, logTeff.p = 4.34±0.02, 7?*,s = 4 .O±O. 57 ? 0 , M.,s = 
7 . 3 ±O. 9 M 0 , and logTeff.s = 4.31 ±0.03 (U05), and Bj.p = 
0.82 kG and Bd,s = 0.5 kG (the minimum values compatible 
with (B^)), we find 77 * > 250 for both stars. Since 77 * >> 
1 , their winds are magnetically confined above the stellar 
surface. The physical extent of the ir magnetospheres, gi ven 
by the Alfven radius Ra (Eqn. lO. lud-Doula et al]l200d i. is 


7?a,p > 4.477?* and Ra,s > 4.267?*. 

In a close binary system it is of interest to compare 
Ra to the physical separation between the stars. U05 found 
semi-minor and semi-major axes of ai sini ~ 0.022 AU and 
02 sini ~ 0.026 AU, or about 0.061 AU and 0.073 AU in the 
orbital plane, assuming an orbital inclination of 21°. This 
is only about 3 7?*,p or 4 7?,,s. Fig. [3] shows a schematic of 
the orbit, in the rest frame of the primary, overlaid with the 
Alfven radii Ra of both stars. We use Alfven radii calcu¬ 
lated upon the assumption of minimum Ba (rather than the 
slightly higher values adopted in Fig. II). There is substantial 
overlap between the two magnetic confinement regions at all 
points during the orbit. Since Fig. [3] involves lower bounds 
for Bd and Ra, this suggests that their magnetospheres 
are very likely to be continuously interacting. If the stars’ 
rotational periods are synchronized, magnetically-ind uced 
spindown JWeber fc Davisiri967l : lud-Doula et al.ll2009l ') may 
well be translated into loss of orbital an gular momentum 
and consequent shrinking of the orbit ([Barker fc Qgilvia 
I 2 OO 9 I I. If the stars are not rotating synchronously, there may 
be magnetic reconnection phenomen a due to the relativ e 
motion of the magnetospheres (e.g. [Gregory et al.l [2014ll . 
Three-dimensional magnetohydrodynamic simulations (e.g. 


[ud-Doula et al.ll2013h will be necessary in order to evaluate 
the importance of such effects. 

There is no evidence o f emission in optical or ultravi¬ 
olet lines ([Petit et al.l 120131 ') . The absence of Hq emission 
is expected, given the weak magnetic fields and winds, and 
the probable slow rotation as inferred from the low values 
of 77 sin 7. The magnetic B2 star HD 3360, which has similar 
magnetospheric properties to those of the e Lupi compo¬ 
nents, shows per iodic UV variability dOskinova et al.[|201l[ : 
[Petit et al.[[2013l l. However, as there is only 1 high-dispersion 
lUE spectrum of e Lupi, UV variability cannot be evalu¬ 
ated. The system is overluminous in X-rays, as expected for 
a magnetic massive star, and indeed is precisely on the X- 
ray lu minosity trend predicted by the semi-an alytic XADM 
model dud-Doula et al.[|2014l : [Naze et al.l[2014h . albeit some¬ 
what harde r than other mag netic, massive stars with similar 
properties dNaze et al.[[2014 1. 

Since the separation between the two stars is relatively 
small, it may be reasonable to consider that the magnetic 
field detected in the secondary is induced by the field of the 
primary. However, the relatively weak magnetic field of the 
primary argues against this: at the distance of the secondary, 
the maximum magnetic field strength of the primary’s dipole 
field is only a few G, assuming a 1/r® decay. This is an 
order of magnitude weaker than the measured (B^) of the 
secondary. In order to produce (B2)~100 G at the distance 
of the secondary, the primary’s magnetic dipole would need 
to be ~15 kG, which is clearly incompatible with the data. 
Moreover, it is unclear how an induced field would explain 
the apparently anti-aligned fields of the two stars. 

The discovery of a doubly-magnetic early-type close 
binary system may generate new insights into our under¬ 
standing of the origin of fossil magnetic fields. The MiMeS 
survey showed MMSs to be relatively rare, approximately 
10% of the populatio n of bright OB stars in the Galaxy 
jGrunhut et al.[[2012h . In remarkable contrast, out of 151 
close binary systems observed in the context of the BinaM- 
IcS survey, not one new MMS has been detected, and in the 
broader sample only one magnetic field (in an F5 star) has 
been detected. This implies an incidence of magnetic stars 
in clo se binary systems of no more than ~2% dAlecian et al.l 

[2OI5II . 

Two scenarios have been prop osed for the origins o f fos¬ 
sil magnetic fields: binary mergers dFerrario et al.l2009li . and 
magnetic flux preserved from the star formation proc ess and 
ampl ified during early pre-main sequence evolution dMestell 
[l999li . 


The merger scenario proposes that the magnetic flux is 
generated by powerful dynamos excited during the merging 
event. The fraction of massive stars expected to undergo 
mergers while on the main sequence is roughly similar t o 
the observed fraction of magnetic stars dSana et al.l [2012 1. 
This scenario seems unlikely to account for the existence 
of a doubly-magnetic massive binary. Two separate mergers 
would be required, but the loss of orbital angular momen¬ 
tum required for one binary pair to merge is likely to be 
transferred to one or both of the other stars, thus widening 
or even destroying the orbit. 

In the second scenario, strong magnetic fields threading 
molecular clouds provide seed fields which are then ampli¬ 
fied during the pre-main sequence convective phase. In this 
case, the historical failures to detect a doubly-magnetic 
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hot star binary are a statistical artifact arising from the 
inherent rarity of MMSs: only ~10% of massive stars are 
magnetic, so only ~1% of close massive binaries should 
contain two magnetic stars. This scenario is supported by 
the similar incidence of magnetic fields, and the similar 
magnetic properties, of magnetic HeAe/Be stars and 
main sequence MMSs, indicating that the presence of 
a fossil field has a lready been e stabli shed during the 
pre-main sequence llAlecian et al.l l2013lf . MHD simula¬ 
tions have also shown that strong m agnetic fields i nhibit 
fragmentation of mole cular clouds JPrice fc Batel l2007l : 
ICommercon et aDl2011^ . Since close binaries are thought 
to be p rimordial, rather than the result of gravitational 
capture (iBonnell fc Batelll994h . the rarity of close magnetic 
binaries may also be explained by this inhibition. The 
discovery of a close binary with two MMSs might be 
thought to be inconsistent with magnetic inhibition of 
fragmentation; however, the relatively weak magnetic fields 
of the e Lupi system may be compatible with an origin in 
a weakly magnetized core, within w hich fragmentation may 
still occur (ICommercon et aLll201ll ). However, this scenario 
is not without difficulties, as simulations of magnetized 
core collapse have tended to show that magnetic flux is not 
preserved (Masson et ah, priv. comm.) 
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